figure;
d_z = data_z(icell, :);
subplot(2,1,1);
hold on;
plotHz(5, d_z);
xlim([0 330]);
barpatch(events{1}(:,1), events{1}(:,2)-events{1}(:,1), [-3 3]);
xlabel('Time (sec)')
ylabel('zscore');
title('to pups');

subplot(2,1,2);
hold on;
y = d_z(zerostate);
x0 = d_z(behavioralstate);
plotHz(5, y);
a = plotHz(5,x0, 'r');
ylim([-3 3]);
xlabel('Time (sec)');
ylabel('zscore');
title(sprintf('ROC=%.3f, p=%.3f', aucs(icell, istate), pvals(icell, istate)));

[h, p] = ttest2(y, x0)
[p, h] = ranksum(y, x0)
